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ABSTRACT 

The hydrogen-deficient [WCL] type central star HD167362 and its planetary nebula 
(PN) SwSt 1 are investigated. The central star has a carbon-rich emission line spectrum 
and yet the nebula exhibits a 10-fim emission feature from warm silicate dust, perhaps 
indicating a recent origin for the carbon-rich stellar spectrum. Its stellar and nebular 
properties might therefore provide further understanding as to the origin of the [WCL] 
central star class. 

The central star optical and UV spectra are modelled with state of the art non- 
LTE codes for expanding atmospheres, from which the stellar parameters are de- 
termined. Using the Sobolev approximation code ISA- Wind, we find Toff=40000 K, 
log(M/M0 yr-i)=-6.72, L=8900 L© (for a distance of 2.0 kpc), and Voo ^ 900 km s'^. 
The abundance mass fractions for helium, carbon and oxygen is determined to be 37%, 
51% and 12%, respectively. From this we derive C/0=4.3 (by mass), confirming that 
the star suffered efficient third dredge-up. The nitrogen abundance is close to zero, 
while an upper limit of <10% by mass is established for H. The model uses a com- 
posite beta velocity law which allows us to reproduce the optical line profiles. The 
overall shape of the de-reddened spectrum agrees with the V-scaled (mv=11.48 mag, 
E(B-V)=0.46 mag) model atmosphere, showing the nebular-derived reddening to be 
consistent with the reddening indicated by the stellar analysis. We confirm our model 
results by using the co-moving frame code CMFGEN, although a few differences remain. 

The PN has a high electron density (log(Ne/cm~'^)=4.5) and a small ionized 
radius (0.65 arcsec - measured from the HST-WF/PC H/3 images), indicating a young 
object. Its nebular abundances are not peculiar. The nebular C/0 ratio is close to 
solar, confirming the PN as an O-rich nebula. The nebular N/0 ratio of 0.08 is not 
indicative of a Type-I PN although the high stellar luminosity points to a relatively 
stellar mass. Near-IR spectroscopy is presented and fitted together with IRAS fluxes 
by using two blackbody curves with temperatures 1200 K and 230 K, indicating the 
presence of hot dust. We also report the first detection of II2 in this young and compact 
PN. 

All of the published spectroscopy since the discovery of SwSt 1 in 1895 has been 
re-examined and it is concluded that no clear spectral variability is seen, in contrast 
to claims in some previously-published studies. If an event occurred that has turned 
it into a hydrogen-deficient central star, it did not happen in the last 100 years. 

Key words: stars: individual: HD167362 - stars: Wolf-Rayet - stars: abundances - 
stars: AGB and post-AGB - stars: atmospheres - stars: evolution - planetary nebulae: 
individual: SwSt 1 - planetary nebulae: general 



1 INTRODUCTION 

[WCL]0 central stars of Planetary Nebulae (PN) are 
hydrogen-deficient central stars that exhibit emission line 
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* 'WC stands for Wolf-Rayet of the carbon sequence, 'L' stands 
for 'late' as opposed to the hotter members of this class, the 
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spectra very similar to massive Wolf-Rayet stars. Until re- 
cently it has been impossible to calculate stellar evolution- 
ary models where the star eliminates all of its hydrogen- 
rich envelop at the top of the Asymptotic Giant Branch 
(AGB). This prompted the association of these objects with 
the born-again scenario (Iben et al. 1983), whereby a white 
dwarf experiences a last helium shell flash, after which it 
rejoins the tip of the AGB for a second time and repeats its 
evolution, this time as a hydrogen- deficient central star. 

The PN around [WC] central stars do not appear to 
have any characteristics which distinguish them from the PN 
of hydrogen-rich central stars. A marginally higher nebular 
C/H ratio (De Marco & Barlow 2001), as well as a slightly 
higher PN expansion velocity (Gorny & Stasinska 1995) , do 
not appear to be sufficient to characterise the sample's evo- 
lution. Other puzzling peculiarities concern spectral appear- 
ance difi'erences amongst members of the hydrogen-deficient 
central star group. Some hydrogen deficient central stars ap- 
pear to have weaker lines than others (called weak emission 
lines stars (wels) by Tylenda, Acker & Stenholm (1993)), 
while in the strong-lined [WC] group, there appears to be a 
classification gap, with [WC] stars populating the [WCll-8] 
class (called [WCL]) and [WC5-W01] (called [WCE]; Pot- 
tasch 1996, Crowther, De Marco & Barlow 1998). Addition- 
ally, while it is currently believed that [WCL] stars evolve 
into [WCE] stars, the stellar abundances seem to differ in the 
two groups, with the more evolved group having a lower car- 
bon abundance than the less evolved group (see De Marco & 
Barlow (2001) for a summary of results from the literature), 
a fact that cannot be easily explained. 

On the other hand. Infrared Space Observatory (ISO) 
spectroscopy of PN with [WCL] central stars has shown 
that all have crystalline silicate features in their mid- to 
far-infrared spectra, in addition to the previously known 
carbon-rich PAH-type features that dominate their mid- 
infrared spectra, betraying the simultaneous presence of 
both carbon-rich and oxygen-rich dust (see Barlow 1997, 
Waters et al. 1998, Cohen et al. 1999, Cohen 2001). The rel- 
atively warm temperatures deduced for the 0-rich dust, to- 
gether with the high nebular densities of their compact neb- 
ulae, appear to imply that the C-rich central stars evolved 
relatively recently from an O-rich composition. This would 
in turn suggest that these [WCL] PN have been on the AGB 
recently and so did not result from a born-again evolution- 
ary phase, for which low-density nebulae are expected. This 
would also imply that the carbon-rich [WCL] stars system- 
atically result from oxygen-rich AGB giants. 

Recent improvements in theoretical stellar evolutionary 
calculations have succeeded in eliminating hydrogen from 
the stellar atmosphere at the tip of the AGB, so that the re- 
quirement that the central star undergo a born-again event 
in order to get rid of its hydrogen is no longer in place. Ac- 
cording to Herwig (2000) a thermal pulse just before the star 
departs from the AGB can eliminate most of the hydrogen, 
while a well-timed post- AGB pulse (probabilistically rarer) 
could account for the extreme hydrogen deficiency observed 
in [WCL] stars. The combination of the ISO observations 



[WCE]. The brackets were introduced by van der Hucht et al. 
(1981) to differentiate Wolf-Rayet central stars from massive cen- 
tral stars. 



and the new evolutionary calculations therefore contributes 
to the revision of theories about the origin of [WCL] central 
stars and for the first time we may have a viable scenario 
that does not require the star to go through a born-again 
episode. 

The PN SwSt 1 around the central star HD 167362 (we 
will give both star and nebula the name SwSt 1) exhibits 
a broad lO^m silicate emission feature (Aitken & Roche 
1982) and has been suggested to be amongst the youngest of 
[WCL] nebulae, possibly having experienced the transition 
from the AGB in the last century. The spectrum of SwSt 1 
was first reported by Fleming (1895), but the first complete 
spectroscopic analysis was carried out by Swings and Struve 
(1940, 1943). The traditional WR C iv A5806/C ill A5696 
line diagnostic ratio positions it within the [WC9] spectral 
class (Crowther et al. 1998 - see also previous classifications 
by Cohen (1975) and by Carlson and Henize (1979)), but 
the weakness of its lines compared to other [WC9] stars 
(e.g. BD-)-303639, He 2-99) or even cooler [WCIO] stars (e.g. 
CPD-56°8032, He 2-113, M 4-18) prompted Crowther et al. 
(1998) to flag it as pecuUar and has therefore essentially left 
it in a class of its own. Tylenda et al. (1993) listed it amongst 
their wels because of the weakness of its spectral lines, but 
SwSt 1 does not resemble stars in this class either (e.g. Cn 3- 
1, De Marco 2000) since it has many more stellar emission 
features than any wels, indicative of a denser wind. 

In this paper we tried to determine whether the physical 
characteristics of the central star and the PN can shed light 
on the reason for these differences, i.e. whether SwSt 1 is at 
an intermediate phase between stars in the [WCL] class and 
their ancestors. In the Sections that follow we present opti- 
cal echelle spectra, near-infrared spectra and archival lUE 
spectra, as well as HST images of this central star and its 
PN. Using these data we carry out a thorough analysis, to 
try and single out the peculiarities of this object. In Sec- 
tion 2 we present our observational data, while in Section 3 
we derive some observational parameters. The distance, lu- 
minosity and mass are discussed in Section 4. The stellar 
analysis is tackled in Section 5, while the nebular analysis is 
carried out in Sections 6 to 8. Finally (Section 9) a review of 
the secular variability is reported and conclusions are drawn 
in Section 10. 



2 OBSERVATIONS 

SwSt 1 was observed with the 3.9 m Anglo- Australian Tele- 
scope (AAT) on May 14 and 15, 1993, using the 31.6 1/mm 
grating of the UCL Echelle Spectrograph (UCLES), with a 
1024x1024 pixel Tektronix CCD as detector. Four settings 
of the CCD were needed in order to obtain spectral cover- 
age in the range 3600-9700 A. At each wavelength setting, 
5 arcsec wide-slit spectra were obtained for absolute spectro- 
photometry, while 1.5 arcsec narrow-slit spectra (R=50,000) 
were obtained for maximum resolution in the range 4100- 
5400 A. The continuum signal-to-noise ratio of the averaged 
exposures ranges from ~5 in the far blue to ~30 in the red. 
A log of the observations is presented in Table 0. 

The data were reduced using the iraf (V2.10)p] pack- 

t IRAF is written and supported by the National Opti- 
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Table 1. Log of the 1993 UCLES observations of SwSt 1. Only 
the second night was photometric. 
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We can also report on the first detection of H2, at 
2.1218 fim (u=l-0, S(l)) from tiiis young, late WC-type 
central star PN system. Although the presence of H2 is well 
established in slightly more evolved [WC] systems such as 
the [WC9] PN BD+30°3639 (Kastner et al. 1996), it has 
not been detected from some young, high-density PN with 
bright central stars, e.g. IC 4997 and IC 418 (Zuckerman & 
Gatley 1988, Kastner et al. 1996). This appears to be due 
to their strong /i'2-dissociating UV radiation fields. The rel- 
ative weakness of SwSt I's H2 emission may indicate that 
the neutral regions surrounding this nebula are also mainly 
atomic rather than molecular. This is in agreement with the 
fact that no CO is detected in the neutral envelope of SwSt 1 
(Huggins et al. 1996). 



age at the UCL Starlink node. Wavelength calibration was 
with respect to comparison Th-Ar arc exposures. Absolute 
flux calibration was with respect to the BSiv star HD 60753, 
whose energy distribution has been measured by Oke (1992). 
Details of the laborious echelle data flux calibration proce- 
dure can be found in De Marco, Barlow & Storey (1997; 
hereafter DBS97). Subsequent data reduction was carried 
out with the Starlink dipso package (Howarth et al. 1995). 
After flux calibration, the narrow-slit spectra were scaled 
up to the continuum level of the wide-slit spectra. 

International Ultraviolet Explorer (lUE) observations 
of SwSt 1 were also used. The retrieved LWR and SWP 
spectra were divided up according to the aperture used for 
their acquisition and the weighted average used. One SWP 
observation (SWP17068) was taken in high resolution mode 
and although its signal-to-noise ratio was rather poor, it 
could be used for parts of our analysis. 

2.1 Near-IR Spectroscopy: First Detection of H2 
in a Young [WCL]-Type Central Star 

Long slit, near-IR spectroscopy of SwSt 1 was acquired on 
Sept 1, 1999 with the NTT, using the Son OF Isaac (SOFI) 
instrument, a 1024x1024 pixel NICMOS detector, and low 
resolution IJ (GRB) and HK (GRR) gratings. The spectral 
coverage was 0.94-1.65 fim and 1.50-2.54 pm, respectively, 
with dispersions of 7.0 A/pix and 10.2 A/pix. The 0.6 arc- 
sec slit provided a 2-pixel spectral resolution of 14-20 A. 
The total integration time was 480 sec at each grating set- 
ting. The removal of atmospheric features was achieved by 
observing HD 169101 (AOV) immediately before and after 
SwSt 1, at a close airmass (within 0.03). Similar observations 
of HD 166733 (F8) permitted a relative flux correction, using 
a T=6100 K model atmosphere normalized to V=9.59 mag. 

A standard extraction and wavelength calibration was 
carried out with iraf, while FIGARO (Shortridge, Mey- 
erdierks & Currie, 1999) and dipso (Howarth et al. 1995) 
were used for the atmospheric and flux calibration, first ar- 
tificially removing stellar hydrogen features from the B3V 
spectrum. Our relatively fiuxed dataset was adjusted to 
match previously published near-IR photometry (Allen & 
Glass 1974) via convolution with the appropriate filter pro- 
files. The spectrum is presented in Appendix ^. 

cal Astronomy Obs ervatories (NOAO) in Tucson, Arizona; 
http: / /iraf.noao.edu/l 
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3 OBSERVATIONAL QUANTITIES 

3.1 Apparent Magnitudes 

We convolved our observed central star wide slit spectro- 
photometry of SwSt 1 with broad-band V and B filter pro- 
files (kindly made available to us by Dr. J.R. Deacon) to ob- 
tain an estimate of the apparent V and B magnitudes of the 
stars. We derive V=11.48 mag and B=11.54 mag. Convolv- 
ing with narrow-band (Smith 1968) filter profiles, we derive 
■i;=11.69 mag and 6=11.73 mag. The hipparcos and ty- 
CHO catalogues quote a Johnson V magnitude of 10.94 mag, 
Hpmag = 10.86 mag (median magnitude in hipparcos sys- 
tem) and an observed Hp range of 10.79-10.97 mag. How- 
ever, SwSt 1 is near the faint limit of the tycho detector 
of ~12 mag. The visual photometry of A. Jones, reveals a 
mean magnitude of (10.4 ± 0.5) mag with no sign of vari- 
ability (Jones et al. 1999). 

We finally decided to adopt our own spectro- 
photometric values since our spectrophotometry aligns per- 
fectly with lUE and near-IR SOFI data and since broadband 
photometry of emission line central stars can be contami- 
nated by stellar and nebular emission line fiuxes. Although 
the visual photometry of Jones et al. (1999) shows the lack 
of intrinsic stellar variability, photometry of SwSt 1 carried 
out by A. Landolt (priv. comm.) on February 19, 20 and 
23, 2001, indicates V=11.57 mag (in the standard system of 
Landolt 1992) ^, slightly fainter than our own spectropho- 
tometry as well as all the other estimates and indicating that 
variability might be present. 

3.2 Nebular Radial and Expansion Velocities 

Nebular radial velocities were obtained by determining the 
average wavelength shift of the nebular Balmer lines using 
a single Gaussian fit. Even when a single Gaussian did not 
fit the data very well, it was felt that this would be the best 
way to locate the overall centre of the line without biases 

t This is the average of six estimates, two per night in the range 
11.55—11.59 mag. The aperture used was 14 arcsec wide. Two 
or three faint stars were included in the aperture, such that this 
estimate might be slightly too high, yet it is fainter than anything 
reported before. From the same photometry the following colours 
are reported: B-V = -1-0.078, U-B = -0.960, V-R = -1-1.283, and 
R-I = -0.498. 
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Table 2. Radial and expansion velocities (Columns 1 and 2) for 
SwSt 1 obtained from the Balmer emission lines. Radial velocities 
for the Na I D absorption line components are also presented 
(Column 3). Finally, distances derived from two different methods 
are presented in Columns 4 and 5. 



LSR RV 



PN Expansion 
Velocity 
(km s~^) (km s"-'-) 

-8.1 ± 1.0 21 ± 2 



LSR RV of 

Na I 
(km s~^) 
-78,-38,+8 

"Distance derived from Na I D positive line components 
'Distance derived from the IR flux method 



D" 



(kpc) 
4.9 



(kpc) 
3.0 



Table 3. Wide-slit (5 arcsec) fluxes (ergs cm~-^ measured 
for the nebular hydrogen lines of SwSt 1, compared with values 
from the literature. 



Reference 
AAT (1993) 
Perek (1971) 
FGC84 
dFV87 
AST89 



H/3 
3.19x10"" 
1.66x10"" 
2.14x10"" 
1.82x10"" 
5.89x10 



H7 
1.20x10"" 

7.24x10"^^ 



H<5 

6.53xl0"i2 



11 



FGC84 = Flower, Goharji and Cohen (1984) 
dFV87 = de Freitas Pacheco and Veliz (1987) 
AST89 = Acker, Stenholm, Tylenda (1989) 



introduced by asymmetry. The heliocentric radial velocity 
was measured to be (-17.8 ± 1.0) km s"^ from an aver- 
age of the Hp, H7 and HiJ lines; the LSR radial velocity 
was determined to be (-8.1 ± 1.0) km s"^. These estimates 
compare well with those of Schneider et al. (1983; Vlsr=- 
8.7 km s"^, VHciio="18.6 km s~^). Radial velocities are listed 
in Table |2|, together with nebular expansion velocities mea- 
sured from averaging the HWHM of HP, H7 and H5 (Hq 
was saturated). 



3.3 Reddening 

The reddening was derived by comparing the relative fluxes 
in the nebular H/3, H7 and H5 lines (Ha was saturated). Us- 
ing the Case B hydrogen recombination coefficients of Storey 
and Hummer (1995) for a nebula with electron density Ne 
= 10* cm"^ and Te = 10* K and the Galactic reddening 
law of Howarth (1983), values of E(B-V) were derived. The 
reddening was also derived from a comparison of the radio 
and H/3 fluxes, using the method of Milne and Aller (1975). 
In Table ^ we list Balmer line fluxes measured from wide-slit 
AAT spectroscopy, together with fluxes from the literature. 
Considerable discrepancies exist between the H/3 flux mea- 
sured here and measurements found in the literature. From 
Table ^ we see that our own value is about one third higher 

Table 4. Reddenings derived from Balmer lines ratios and from 
the radio-H/3 baseline for SwSt 1. 

Ratio E(B-V) E(B-V) E(B-V) E(B-V) 

AAT (1993) FGC84 dFP87 P71 

Ha/H/3 - - - - 

H/3/H7 0.51" - - - 

H/3/H<5 0.41" - - - 

radio/H/3 0.30 0.34 0.41 0.50* 

FGC84 = Flower, Goharji and Cohen 1984 
GD85 = Goodrich and Dahari 1985 
dFV87 = de Freitas Pacheco and Veliz 1987 
P71 = Perek 1971 

" We adopted 0.46, the average of 0.51 and 0.41 
^ Derived by us using his H/3 flux 



than most of the other values, although Acker, Stenholm and 
Tylenda (1989) measured an even higher H/3 flux. The pos- 
sibility of a true variation in time was therefore investigated, 
but no clear trend exists in either the H/3 fluxes or in the 
radio flux measurements. Although these large differences 
remain suspicious, we must attribute them to calibration 
and measurement errors and possibly to some slit/aperture 
widths not including all of the PN. 

In Table ^ we list our results together with reddenings 
selected from the literature for comparison. From the H/3- 
H7 baseline in our UCLES wide-slit (5 arcsec) spectrum, we 
determine E(B-V)=0.51, while from the H/3-H(5 baseline we 
find E(B-V)=0.41. Radio fluxes at 15 GHz have been mea- 
sured by Kwok et al. (1981; (207 ± 11) mjy), Milne and Aller 
(1982; (240 ± 12) mJy) and Aaquist and Kwok (1990; (175 
± 17) mJy). The weighted average ((211 ± 17) mJy) corre- 
sponds to a 5 GHz flux of 243.0 mJy for optically thin free- 
free emission, which we adopted, together with our UCLES 
value for the H/3 flux, to derive E(B-V)=0.30. A similar 
value is also obtained by nuUying the 2200 A feature. Com- 
parison of the Perek (1971) photo-electric flux with the mean 
radio flux yields E(B-V)=0.50. 

We adopt E(B-V)=0.46±0.05, the mean obtained from 
the H/3/H7 and H/3/H(5 ratios, as the most self-consistent 
estimate we have. Pre-emptying the results of Sec. the 
slope of the synthetic stellar atmosphere model carried out 
with the ISA- Wind code, which flts the stellar emission line 
intensities, agrees with the nebular line reddening estimate 
of E(B-V)=0.46. For the model carried out with the cmf- 
GEN code, the best match is obtained with E(B-V)=0.48. 
These values are also in reasonable agreement with the de- 
termination of E(B-V)=0.50 quoted above, from the radio 
and the photo-electric H/3 fluxes. 

Finally, the near-IR spectrum is dominated by strong 
nebular emission lines of H I (e.g. P/3, Br7) and can be used 
to obtain a further confirmation of the reddening. For E(B- 
V)=0.46 and the nebular properties derived in Sec. |^, agree- 
ment with Case B recombination theory (Storey & Hummer 
1995) is very good. For instance, J(P/3)//(H/3)=0.16 (0.15 
from Case B), or /(Br7)//(H/3)=0.024 (0.024 from Case B). 
Pq is an exception, but this lies in a region of poor atmo- 
spheric transmission. 



4 LUMINOSITY AND DISTANCE 

In Table ^ we summarise the distances that have been de- 
termined for SwSt 1 in t he li terature, together with our two 

and 4.2. Below we explain our 



estimates from Sections 4.1 



two methods but admit that there is no convincing estimate 
for the distance to SwSt 1 beside the fact that after iterat- 
ing between an assumed distance and our models we have 
gained further insight which restricts the range of possible 
distances (see Sees. ^ and P), leading us to finally adopt a 
distance of 2.0 kpc (Section 4.3). 



4.1 Distance from a Magellanic Cloud Luminosity 
Calibration 

DBS97 used Magellanic Cloud Wolf-Rayet central stars to 
derive a mean central star mass for WC central stars of 
0.62 Mq. This mass was then applied together with the 
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Table 5. Summary of the distances determined for SwSt 1 in the Uterature using various methods. 



Distance Method 
(kpc) 

1.0 extinction method 

1.1 modified Shklowskii 
>1.2 HIPPARCOS 

1.4 Shklowskii 

1.4 M„ - Tt relationship 

3.0 LMC [WC] star calibration 

3.1 M and gravity 

3.2 I(H;9) for opt. thick PN" 

3.8 ionized mass - R(PN) and radio Tb - R{PN) corr 
4.7 radio cont. brightness and R(PN) 

4.9 Galactic rotation curve 
5.0 L* 

''We used their Eq. 11 and our log{I{H/3))=-10.1 from Sec. 



Reference 

de Freitas Pacheco & Veliz 1987 

de Freitas Pacheco & Veliz 1987 

Acker et al. 1998 

Cahn et al. 1992 

de Freitas Pacheco & Veliz 1987 

this work 

Zhang 1993 

Kingsburgh &; Barlow 1992 
Zhang 1995 

van de Steene and Zijlstra 1994 
this work 
Zhang 1993 



helium-burning tracks of Vassiliadis & Wood (1994) to de- 
rive luminosities for CPD-56°8032 and He 2-113 corre- 
sponding to their effective temperatures, which were in turn 
used, together with their integrated IR fluxes, to determine 
their distances. The PN of SwSt 1, with a high electron den- 
sity, a small apparent radius and a significant IR excess ap- 
pears to be a good candidate to apply the same method. We 
have obtained from the ISO archive the SWSOl spectrum of 
SwSt 1 (from the programme of S.K. Corny, TDT47101511, 
reduced with SWS OLP Version 9.5, Szczerba et al. 2001). 
This spectrum is plotted in Fig. ^ together with its pho- 
tometric fluxes at 1.65 and 2.2 /im (Allen & Glass 1974) 
and colour-corrected IRAS fluxes. Integration under the IR 
spectrum of SwSt 1 between 1.6 and 100 fim results in a flux 
of 1.48x 10~* erg s~^ cm~^, corresponding to an infrared 
bolometric luminosity of 463 D^(kpc) Lq. For Teff ~40000 K 
(Section the helium-burning tracks of Vassiliadis and 
Wood (1994) predict a luminosity of 4280 Lq, where we 
interpolated between the tracks for 0.600 Mq and 0.634 
Mq cores. The distance obtained in this way is 3.0 kpc (a 
Teff~35000 K predicts 4670 Lq and a distance of 3.2 kpc). 

The three key assumptions of this method are: (a) 
SwSt 1 has the same mass as the LMC [WC] stars, (b) all 
its luminosity is re-radiated in the IR by dust and (c) the 
helium-burning track calculations are valid. By far the most 
compromising assumption is that all [WC] stars have the 
same mass (as suggested by De Marco and Crowther (1998) 
this is likely not to be the case). We have therefore deter- 
mined the distances corresponding to all three masses cal- 
culated by Vassiliadis and Wood (1994), namely 0.57, 0.60 
and 0.63 M0. These are 2.1, 2.5 and 3.6 kpc, respectively. 
We therefore conclude that if assumptions (b) and (c) hold, 
the distance is in the range ~2.0-3.5 kpc. 



4.2 Galactic Rotation Curve Distance 

In order to obtain an alternative value for the distance, 
we applied the insterstellar sodium D line method used by 
DBS97. In Fig. ^(a) we show the spectral range containing 
the Na I D lines, while in Fig. |l](b) we show Na I Dl and 
D2 lines in their respective LSR rest frames. From Fig. |l|(b) 
and (c) we see that the line of sight to SwSt 1 has a positive 
rotation curve, while most of the Na I D line components 
have negative radial shifts. The only interstellar component 



withpositive velocity is measured at (+8 ± 2) km s~^ (see 
Fig. p|(b)) and corresponds to a distance of (4.9 ± 0.6) kpc, 
where the uncertainty derives solely from locating the cen- 
tre of the trough. The underlying stellar emission lines (C ll 
AA5889,92, C lll A5894 and He ll A5897) are too weak to 
achieve a reliable rectification, which would lead to an in- 
crease in the measurement precision. Moreover the Na I D 
lines have emission components (as noticed by Dinerstein 
et al. (1995)) and a multiple cloud fit would over-interpret 
the data, so we limited ourselves to measuring the centres 
of the absorption components. The negative radial velocity 
troughs are likely to derive from the PN's neutral envelope 
and are discussed in Appendix ^. 

4.3 The Adopted Distance 

Although each method listed in Table ^ has its limitations 
and it is beyond the scope of this paper to review each 
of the different estimates, we should note that two of the 
best methods to estimate distances to Calactic objects suf- 
fer from additional caveats. SwSt 1 is too far for the HIP- 
PARCOS distance, derived from a parallax of (8.9 ± 6.0) mas 
(HIP89535; Acker et al. 1998), to be considered reUable. Sec- 
ond, the extinction distance is based on the assumption that 
within 2 kpc of the sun, light suffers 0.4 mag of extinction 
per kpc, and is not based on an extinction map, which would 
make the method more reliable. 

Although our estimates seem to point to a fairly large 
distance (i.e. ~3-5 kpc), we found in our modelling of the 
stellar spectrum that a very high luminosity and central star 
mass were implied by a distance larger than 2 kpc and would 
have required the central star to have evolved significantly in 
effective temperature in the last 100 years. Such evolution 
would have been mirrored by significant spectral changes. 
This however has not been observed (see Sec. We there- 
fore adopt 2 kpc as the distance to SwSt 1, but admit a 
substantial uncertainty. 



5 STELLAR ANALYSIS 

In this Section we carry out a stellar atmosphere analysis 
of SwSt 1 to derive the stellar effective temperature and 
mass-loss rate as well as its carbon and oxygen abundances. 
The line-blanketed stellar flux resulting from this analysis is 
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then used as input for the photo-ionization model of the PN 
(Sec. |). 

5.1 The Near-IR, Optical and UV Stellar Spectra 

The spectrum of SwSt 1 contains a mixture of stellar and 
nebular emission lines. The stellar lines are weaker and fewer 
than for the [WCIO] central stars CPD-56 80° 8032 and He 2- 
113 (DBS97), while the nebular lines are stronger. The neb- 
ular Balmer lines, the lines of He I, and, to a lesser extent, 
those of C II and O ll, could be blends of stellar and nebular 
emission. In particular in the UV C ill] A1909 is a blend 
of stellar and nebular emission, with the nebular emission 
contributing about 40% of the flux to this feature (Fig. |^. 
The C ll] A2326 feature is thought to be mostly of nebular 
origin, although this cannot be checked since only low res- 
olution lUE data exists for this wavelength range. The UV 
C lll-iv and Si iv e mission lines are entirely of stellar origin. 

In Sect. [3.I.1I , we will show that it is unlikely that hy- 
drogen is present in the star, so that the main contributors 
to the Balmer lines are likely to be nebular hydrogen and, 
to a much lesser extent, stellar He 11. Additionally, there is 
no evidence of a significant stellar contribution to the strong 
nebular He I lines. Except for A5876, which exhibits a broad 
pedestal at the base of the nebular component, other strong 
He I lines at 3889, 4471, 4713 and 6678 A, show no evi- 
dence of having a stellar component. The broad pedestal at 
the base of He l A5876 is likely to belong in part to C ill 
M 20, with components at (in order of intensity) 5894.07, 
5871.69, 5858.35, 5880.54, 5863.24 and 5872.10 A, since the 
four strongest components are observed at approximately 
the predicted relative strengths. In the near-IR we observe 
He I at 1.083 fira, although the lower resolution of this spec- 
trum does not allow us to distinguish between nebular and 
stellar components to this line. Stellar He 11 A4686 is clearly 
visible in emission. He 11 A5411 has a weak P-Cygni absorp- 
tion but its emission, if present, is blended with the nebular 
[Fe III] line at 5412.0 A; He 11 A4541 has a weak emission to 
a blue-shifted absorption. 

C III lines dominate over C 11 lines in SwSt I's spectrum, 
indicating a higher degree of ionization than for the two 




-2000 -1000 -500 500 1000 1500 2000 2500 

Velocity (km/ s) 

Figure 2. The C IV resonance doublet proflle in the lUE high res- 
olution spectrum (SWP17068HL). The x-axis has been converted 
to velocity space; v=0 km/s corresponds to 1548.2 A, the bluer 
of the two components of the doublet (at 1548.2 and 1550.8 A). 
The vertical line marks the chosen vijiack=^900 km . 



[WCIO] stars analyzed by DBS97. In the optical we observe 
C III A4650, A5696, A8500. A8500 is a blend of three lines: 
the stellar C ill A8500, a strong nebular line with radial 
velocity-corrected wavelength at 8502.1 A, likely to be itself 
a blend of H I Paschen 16 and [CI ill] A8500 ([CI ill] A8481 is 
also present), and a weaker line at 8499.7 A. An unidentified 
nebular line at this spectral position was also observed by 
Keyes et al. (1990) in the spectrum of NGC 7027. In the 
near-IR SOFI spectrum we observe C ill at 9710 A, and at 
1.198 ^m. 

C IV AA5801,12, A5470 and A4441 are observed, al- 
beit weakly, while in the UV we see C iv AA1548,1551. 
C II lines are observed at 4267 A, 6578 A and 7231.3 A 
(blended with telluric absorptions). All C 11 features exhibit 
stellar as well as nebular components. Semi-forbidden C 11] 
AA2324.2,2325.4 (likely to be entirely of nebular origin) and 
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C III] AA1906.7,1908.7 (having a nebular as well as a stellar 
component - Fig. ^ , are also observed. 

O III A5592 is observed in emission with a weak P-Cygni 
absorption, while the O ll stellar spectrum is very weak (e.g. 
AA4591. 0,4596.0). No nitrogen lines are observed, although 
there is a possibility that a feature observed in the lUE high 
resolution spectrum at 1718.92 A is indeed N iv A1718.5 
and not Si iv A1722.5. The optical lines of nitrogen used 
by Leuenhagen & Hamann (1998) (i.e. N ill AA4097,4104) 
to determine an upper limit for the abundance of nitrogen 
are not observed in our spectrum. N ii A3995 and A4682 
are also unobserved. The N v line marked by Flower et al. 
(1984) in the UV spectrum at 1438.8 and 1442.8 A is almost 
certainly C ill A1347.4. Si iv AA1393. 7,1402.8 are observed 
in the UV spectrum, as are weak optical lines at 4089 and 
4116 A (Ml). 

From the high resolution lUE spectrum of the C iv 
A1550 line we measure a Ubiack (which corresponds to the 
wind terminal velocity Voo', Prinja, Barlow & Howarth 1990) 
of (900 ± 150) km s~^ (Fig. ^). We notice how this value is 
much larger than anything implied by the optical line widths 
and P-Cygni features (which is why Leuenhagen & Hamann 
(1998) adopted a much lower value for v.oo of 400 km s~^). 
This is likely to be a consequence of a particularly low wind 
density. 



5.1.1 The B aimer Lines 

DBS97 noted that the Balmer line shapes of the PN around 
the two [WC] central stars CPD-56°8032 and He 2-113 have 
a broad pedestal. Leuenhagen et al. (1996) determined an 
upper limit for the stellar hydrogen abundance by fitting this 
pedestal with synthetic profiles from their non-Local Ther- 
modynamic Equilibrium (non-LTE) code. However DBS97 
pointed out that the same broad pedestal was present at the 
base of forbidden lines which are entirely of nebular origin. 
This indicates that either a complex nebular velocity field 
is the cause of the irregular nebular line shapes and of the 
pedestals at the base of the Balmer lines, or perhaps more 
likely, that the pedestals observed at the base of these strong 
lines is due to the grating instrumental profile. 

Since the presence of hydrogen in cool [WC] central 
stars would give a clue about their origin, we have inves- 
tigated the nebular lines of SwSt 1. In Fig. ^ we present 
Gaussian fits to Hq and H/3, along with other forbidden 
lines. Although the Balmer lines do present a base which is 
broader than the Gaussian fit, many other nebular lines, i.e. 
[N ll] AA6548,6583 and [O lll] A5007, have a similar shape. 
There is therefore no conclusive evidence for the presence of 
stellar hydrogen. In Sect. we demonstrate how a hydro- 
gen abundance of 10% (by mass) as determined by Leuen- 
hagen and Hamann (1998) is not inconsistent with our own 
fits. 



radiation transfer is implemented using an improved ver- 
sion of the Sobolev approximation, which has been shown 
to be appropriate for the thick atmospheres of Wolf-Rayet 
stars (de Koter et al. 1997, Crowther et al. 1999 and De 
Marco et al. 2000). 

The wind of SwSt 1, however, is less dense than the 
winds of most Wolf-Rayet stars, as can be inferred from the 
overall weakness of its spectral lines. This indicates that the 
inner, slow moving parts of the stellar atmosphere are ex- 
posed in this region, and the Sobolev approximation might 
break down. To check on the results obtained with ISA- 
Wind, we used the co-moving frame code CMFGEN (Hillier & 
Miller 1998). CMFGEN confirmed the results obtained with 
ISA- Wind. However, some of the differences between the two 
codes are instructive, and although in what follows we re- 
fer mostly to ISA- Wind, results obtained using CMFGEN are 
presented alongside. 

Line blanketing is included in the ISA- Wind model by 
a Monte-Carlo approximation, which includes millions of 
spectral lines. This is described in detail by Schmutz et al. 
(1991). Comparison tests have been carried out (Crowther 
et al. 1999, De Marco et al. 2000) and show the technique to 
compare reasonably well to a detailed line-blanketing cal- 
culation carried out with CMFGEN (Hilher & Miller 1998), 
although the flux shortwards of 400 A was too hard in the 
Monte-Carlo technique. 



5.2.1 The Model Atom 

The number of atomic levels included in the calculation in- 
fluences the line strengths. Ideally one would want to employ 
as many lines as possible of as many elements as possible. 
Unfortunately the larger the atomic model the slower the 
calculation. As a result a compromise has to be adopted. 
The atomic model used by our ISA- Wind calculation uses 
Opacity Project data (Seaton 1995) and is summarised in 
Table B The use of super-levels (SL) reduces computational 
times while at the same time including high-lying levels. 
Super-levels are groups of levels within which populations 
are calculated in LTE. The grouping scheme changes for dif- 
ferent ions. In Table ^ we indicate the scheme used: "n=4" 
means that all levels up to and including principal quantum 
number 4 are included, "sing/trip" means that singlets and 
triplets are grouped separately; "thr." indicates that enough 
super-levels were constructed to include all levels listed in 
the Opacity Project database up to the ionization thresh- 
old. The effect of slight variations in our grouping scheme 
(as well as the total number of SL) on the synthetic lines 
was studied for the cases of C ll and C ill and found not 
to alter the synthetic lines. Only helium, carbon, oxygen, 
silicon and iron are all included explicitly in the CMFGEN 
calculation, but in somewhat greater detail than in the case 
of ISA- Wind. 



5.2 Description of the Model Codes 

The code adopted to carry out the stellar analysis is the 
improved Sobolev approximation code ISA- Wind (de Koter 
et al. 1993, 1997). This modelling program solves the ra- 
diation transfer equation under the condition of statistical 
equilibrium in a spherical, extended atmosphere. The line 



5.3 Photon Loss 

"Photon loss" is the name given by Schmutz (1997) to the 
interaction between the photon fleld of the He ll Lya line at 
303.78 A and nearby metal lines (e.g. the O ill Bowen lines 
at 305.72 and 303.65 A). This can be responsible for a de- 
crease in the ionization equilibrium of the wind with respect 
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Table 6. Summary of the atomic model adopted in our ISA- Wind 
non-LTE simulations of SwSt I's stellar wind. 
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to the equilibrium obtained when no such line-line interac- 
tion is accounted for. In order to fit the observed spectrum 
when photon loss is included in the model, the model's stel- 
lar temperature has to be higher, which results in a higher 
spectroscopic luminosity. This contributes to a lower wind 



performance number and allowed Schmutz (1997) to calcu- 
late a velocity structure for the outer wind of WR6. 

De Marco et al. (2000) implemented this approximation 
in the Sobolev code ISA- Wind and showed that results were 
similar to those obtained using the code CMFGEN (for the 
WC8 stellar component to the 7 Velorum binary), where 
line-line interaction is taken into consideration explicitly. 
Here, we have used the "photon-loss" option, which resulted 
in only a slightly higher effective temperature (^2000 K). 
The efi'ect has been shown by Schmutz (1998) not to be 
important in the case of O stars, which have weaker winds. 
This also seems to be the case for the weaker wind of SwSt 1. 
CMFGEN naturally includes line-line interaction and there- 
fore "photon-loss" is automatically taken into account (De 
Marco et al. 2000). 

5.4 Stellar Modelling Technique 

Traditionally, modelling of WR stars starts with matching 
a subset of the observed He I and He 11 lines with model 
line profiles, from which approximate values for the efi'ective 
temperature and mass-loss rate of the star are determined. 
The metal abundances are then varied, to afford fits to a 
representative sample of their lines. Even in modern calcu- 
lations the velocity law is seldom varied from either a simple 
or a composite beta- law (e.g. Hamann & Koesterke 2000). 

A weakness of this technique is that in a significant 
number of cases, the choice of diagnostic lines determines 
what parameters are derived. If excellent fit quality can be 
obtained for a large number of lines, this gives confidence 
in the derived parameters, but if only a small subset of the 
observed lines can be fit by any given model, the param- 
eters are only as good as the choice of diagnostics. Addi- 
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tionally, running non-LTE line-blanketed codes is very time 
consuming, so that once a reasonable match to a set of lines 
is achieved, different parameter spaces are seldom explored 
to determine whether other parameter sets fit the emission 
lines to a similar degree of accuracy. 

The fast calculation speed of the ISA- Wind code allowed 
De Marco et al. (2000) to show that similar synthetic spectra 
can sometimes be obtained with different combinations of 
temperature and mass-loss rate. When only a poor fit quality 
is achieved, or when only a subset of the diagnostic lines is 
fit simultaneously, it might be hard to choose which model 
best represents the star, leading to a degeneracy. 

The spectrum of SwSt 1 presents a formidable challenge 
for spectral synthesis, in that it is impossible to fit all the 
diagnostic lines simultaneously. As a result, it is possible 
to derive different model parameters depending on which 
diagnostic lines we decide to include. The exploration of a 
vast parameter space (stellar temperature, mass-loss rate 
and metal abundance), made possible by our fast code, to- 
gether with the use of observational constraints not usually 
coupled with non-LTE analyses, made it possible to gain a 
reasonable confidence in our final choice of stellar param- 
eters. Our choice of diagnostic lines comprised all lines in- 
cluded within the levels of our model atom, leaving out those 
between high-lying levels which might be under-populated. 
For helium we used He ll A4686 (A5411 is blended with neb- 
ular a [Fe ll] line, and can only be used as an upper limit) 
and He i A5876,4471,6678. As pointed out above, there is 
extremely little He l in this star, if any, so that the lines 
(severely blended with nebular components) could only be 
used to impose an upper limit on the strengths. To constrain 
the ionization equilibrium one has to appeal to carbon: C ll 
A4267, C III A5696 and A8500, and C iv AA5801,12 and A5470 
were employed. These lines are also used to determine the 
carbon abundance. To constrain the oxygen abundance we 
used O III A5592 (all other lines are too weak or dominated 
by the nebula). 

Once the diagnostic lines were chosen, an approximate 
value for the temperature and mass-loss rate was obtained, 
for which reasonable fits were achieved. Since no combina- 
tion of Tt, M and carbon abundance fitted the spectrum 
to a satisfactory level, in that the synthetic lines were al- 
ways too broad in the core and too narrow at the base, we 
experimented with different velocity laws. In particular we 
tested velocity laws with progressively slower flows in the 
inner regions, so as to reduce the line widths, maintaining 
the observationally-derived value of v^a- The electron den- 
sity change induced by the different laws changes the wind 
opacity, with a repercussion on the line ionization equilib- 
rium of the star, and therefore on the line intensities. Con- 
sequently different basic parameters are implied by different 
velocity laws. 

5.5 Stellar Model Results 

Model fits to the stellar spectrum of SwSt 1 are presented in 
Fig. ^ (thick and thin solid lines for ISA- Wind and CMFGEN 
models, respectively, dotted line for the observations), while 
the model parameters are presented in Table 0, together 
with the results of Leuenhagen & Hamann (1998), scaled 
to our adopted distance. In what follows we discuss the re- 
sults obtained with the modelling code ISA- Wind, followed 



by a discussion of additional results obtained with the code 
CMFGEN to later. 

5.5.1 ISA- Wind 

He II A4686 is well reproduced, while He ll A5411 and He I 
AA5876,6678, and other lines which suffer blending with neb- 
ular emission, are within reasonable limits. Other synthetic 
He II lines (e.g. A4542 and A4200) are in absorption with 
wide wings. Since these wings are not observed, and He I 
lines are mostly of nebular origin, we have have not used 
these lines as diagnostics. 

C IV A5801,12 is over-predicted, C iv A5470 is weU 
matched, while C iv A4441 is somewhat weak. In order to 
fit C IV AA5801,12 the effective temperature could not be 
higher than 35 kK, no matter what other parameters are 
adopted. However at this temperature the A5470 line is in 
weak absorption. C iv A5470 has a better-understood for- 
mation mechanism compared to AA5801,12 (Hillier & Miller 
1999, Dessart et al. 2000) and was therefore given priority 
in the fitting procedure. Additionally, at 35 kK the fit to the 
C III A4650 line is much worse, with the synthetic line too 
weak. C III A5696,8500 are also well matched by the model 
but C II A4267 is poorly matched, with the model line being 
much too weak. However the C ll AA6578,6583 blend has 
a very reasonable fit (C ll AA7231,7236 is severely affected 
by telluric absorption and blending with nebular emission 
line components and was therefore not fitted) . The problem 
with C II A4267 might be linked to insufficient levels in the 
atomic models. However we have experimented with more 
levels and included a considerable number of super-levels 
and dielectronic recombination lines. Since a much stronger 
synthetic A4267 (obtained by lowering the effective temper- 
ature) would lead to over-fitting AA6578,83, we still suspect 
the model atom is the reason for the problem. If a cooler 
stellar temperature is used to try and enhance C ll A4267, 
all fits deteriorate, even if the mass-loss rate is reduced to 
restore the He l/He ll equilibrium. Increasing the mass- loss 
rate so as to enhance recombination of C ill to C ll, and so in- 
crease the strength of the C ll A4267 line, also fails. A higher 
mass-loss rate is also excluded by the fact that by increas- 
ing its value by 0.1 dex from the value reported in Table ^ 
all model lines of C ill acquire deep P-Cygni absorptions, 
which are not observed. Additionally the He I lines increase 
in strength, in a way that cannot be compensated for by a 
larger temperature without losing the C lll/C iv balance. 

Increasing the carbon abundance, in steps, from 
C/He=0.2 to 3 (by number), while adjusting other param- 
eters, does not achieve an improvement in the fits. On the 
other hand, once the effective temperature and mass-loss 
rate have been determined, the C/He ratio can be deter- 
mined with a 50% accuracy. Even with the highest carbon 
abundance, the C ll line did not increase for the temperature 
regime which satisfied the other lines. Additionally, when 
enhancing the C/He ratio the ionization equilibrium of the 
star is upset in such a way that He ll lines are drastically 
decreased, which results in having to increase the tempera- 
ture, at which point other carbon lines do not fit any more, 
most noticeably C iv AA5801,12, whose excitation is almost 
entirely dependent on temperature. 

Returning to the C ill line blend at 4650 A, and to line 
shapes in general, we have used a velocity law which is slower 
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Table 7. Stellar parameters from our non-LTE models (iSA-Wind and CMFGEn), compared with those of Leuenhagen 
LH98), before and after scaling them to our own adopted distance. 
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" Limit from model profile fitting does not include possibility of blend with nebular emission (see text). 
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than 13=1 in the inner parts of the wind (in accordance with 
that found by Hillier & Miller (1998), Crowther et al. (1999) 
and De Marco et al. (2000)). For every parameter set, we 
tested a number of velocity laws with progressively slower 
winds in the inner parts of the flow. These are shown in 
Fig. ^. The slower the velocity law, the stronger and nar- 
rower the lines become (for identical stellar parameters) . As 
a result each adopted velocity law implied a new set of pa- 
rameters, resulting in a trend of a somewhat cooler stellar 
temperature and lower mass-loss rate for slower velocities. 
For the slower velocity laws, narrower line profiles were ac- 
companied by more redward-displaced P-Cygni absorptions. 
As a result it was relatively easy to choose the best velocity 
law from the position of the P-Cygni absorption (as well as 
the overall line width and shape) of lines such as C ill A4650 
and C IV AA5801,12. The velocity law which matched the po- 
sition of the P-Cygni profiles is marked as a thin solid line 
in Fig. ^ (left panel) and can be parametrised following the 
formula of Hillier and Miller (1998) with Uoxt = 500 km s"\ 
Voo = 900 km s~^, /3i = l and /32=8. This resulted in the best 
match to the line shapes. In Fig. ^ (right panel) we present 
profile fits obtained using a P=l law, and two slower laws. 
The model parameters were adjusted in each case to repro- 
duce the line's peak intensity. As can be seen, the (3=1 law 
generates profiles which are too broad, while the dot-dash 
line produces profiles too narrow at the base but too broad 
in the core. 

In summary, although the fits presented in Fig. ^ (thick 
solid line) are not excellent, there is no other parameter 
combination of mass-loss rate, temperature and C/He num- 
ber ratio^ which would satisfy the observed spectrum to a 
higher degree of accuracy. The derived quantities and their 
errors are therefore the following: T. = (40 000 ± 2000) K, 
log(M/MQ) = (-6.72 ± 0.1), C/He = (0.46 ± 0.20) and 
O/He = (0.08 ± 0.02) (by number). The wind performance 
number r) (=L/{Mvaoc) is 0.9, indicating that the wind of 
SwSt 1 can be driven radiatively in the single scattering 
limit. The derived stellar luminosity of 8900 Lq is towards 
the high end for a PN central star (implying a relatively large 
stellar mass) . It is difficult to determine the mass which cor- 
responds to this luminosity since the highest core mass for 
which Vassiliadis and Wood (1994) calculated Milky Way 
initial abundance helium-burning tracks was 0.63 M0, but 
inspection of their Fig. 8 indicates that our derived luminos- 
ity for SwSt 1 at 40 000 K is only slightly above the LMC- 
abundance evolutionary track for a 0.679 Mq core mass. 



§ The parameter space explored spanned from 30kK to 45kK in 
stellar temperature, log(M/M0)=-6.O to -6.7 and C/He=0.2 to 
3.0. 



The helium, carbon and oxygen abundances are 37%, 
51% and 12% by mass, respectively, implying a C/0 mass 
ratio of 4.3. The oxygen abundance in particular is obtained 
from fitting the O ill line at 5592 A. The lower and upper 
limits for the O/He ratio allowed by the fit of the A5592 line 
are 9% and 15%. Additional uncertainty derives from fitting 
a single line. We note that if a lower effective temperature 
were to be adopted, the oxygen abundance determined by a 
fit to this O III line would be higher still. 

In Fig. ^ we show the ISA- Wind synthetic spectrum, 
compared with the IR, optical and lUE de-reddened data. 
A slight over-compensation of the 2200 A feature by the 
adopted E(B-V)=0.46 is not considered worrying in view 
of the known variations in the relative strength of this fea- 
ture (Mathis 1994), since the overall continuum slope is well 
matched by our V-scaled synthetic spectrum (if we make an 
exception for the far blue and red edges of the optical spec- 
trum). By adopting reddening values as low as 0.30 (which 
nulls the 2200 A feature) the continuum slope is clearly 
too gentle to match the model, while for a reddening value 
of 0.50 the slope will be too steep. In Fig. ^ we also plot 
the near-IR SOFI spectrum (scaled to the J magnitude of 
Allen & Glass (1974)) and colour-corrected IRAS photom- 
etry values, together with blackbody curves with tempera- 
tures 1200 K and 230 K. These curves indicate that some of 
the dust around this PN is very hot. This is in agreement 
with the study of Zhang & Kwok (1991). 

5.5.2 CMFGEN 

The fit quality obtained with the code CMFGEN is very simi- 
lar to that obtained with ISA- Wind. Although some lines are 
better reproduced (e.g. C iv AA5801,12) others are some- 
what worse (e.g. C iv A4441). The broad He 11 absorption 
lines wings obtained with ISA- Wind are somewhat narrower 
in the models fits of CMFGEN, in better agreement with the 
data. C III A1909 is weaker and more in agreement with the 
observations (note that in Fig. ^ this line is a blend of stellar 
and nebular components, as can be appreciated by inspect- 
ing Fig. ^) . On the other hand the shape of the C lll-iv blend 
at ~4650 A is somewhat worse in the CMFGEN fits, and the 
C III line at 8500 A has an unobserved P-Cygni absorption. 
The fit to the only oxygen diagnostic line, i.e. O ill A5592, 
is similar to the one obtained with ISA- Wind. The stellar 
parameters obtained with cmfgen are presented alongside 
those obtained with ISA- Wind in Table ^ Overall, the de- 
rived stellar parameters for SwStl using CMFGEN match 
closely with those from ISA-wind, except that the stellar 
radius is smaller, implying a higher stellar temperature of 
46kK. The mass- loss listed in Table ^ does not consider 
clumping. Almost identical fits can be obtained with a fiU- 



© 1999 RAS, MNRAS 000, 



12 



0. De Marco et al. 




R(R. -1) Wavelength (A) 

Figure 5. Left panel: alternative velocity laws compared to a /3=1 law (thick solid line). The adopted law (thin solid line) significantly 
improved the line fits, making lines broad at the base and narrow in the core, and reproduced the position of their P-Cygni absorptions, 
while the other laws produced lines which were too broad or too narrow. Right panel: ISA- Wind profile fits to C III A5696 (top) and 
C IV AA5801,12 (bottom). The dotted line represents the data, the thick solid line represents the model calculated using a /3=1 velocity 
law, the thin solid line represents the model with the velocity law used for the final model of this paper, and the dashed line shows the 
slowest velocity law attempted. Model parameters are adjusted to fit the peak intensity of C IV AA5801,12. 



ing factor of 10% with clumping starting aX v — 100 km s~^. 
In this case the mass-loss is 7x10~*Mq yr"'^ or a factor of 
two lower than obtained without clumping. Due to the rel- 
atively low wind density, the electron scattering wings can- 
not be used to constrain the clumping. The only noticeable 
changes are in C ill A9710 and A1909 where clumped models 
give 20% weaker emission. However, since the latter is dom- 
inated by the nebula and the former is only included in our 
low resolution SOFI data, these cannot provide any useful 
constraint. 

The only considerable difference between the modelling 
efforts with ISA- Wind and CMFGEN is that the stellar abun- 
dances are discrepant (Table |^. This is possibly due to the 
different cocktail of atomic data used by the model codes, 
although, as it is frequent with complex non-LTE codes, it 
might also be due to a mix of factors hard to trace. In par- 
ticular the carbon and oxygen mass fractions are 32% and 
15%, respectively, compared to ISA-Wind's 51% and 12%. 
This translates into a C/0 mass ratio of 2.1 for CMFGEN ver- 
sus 4.3 for ISA- Wind. Both CMFGEN and ISA- Wind values are 
higher than solar (~0.4). The maximum value that the C/O 
ratio can attain according to the models (Herwig 1999), is 
the intershell value, ~4. So, although both values are accept- 
able, the IS A- Wind value would imply an extremely effective 
third dredgeup. 

We obtain a good fit of the CMFGEN model flux distribu- 
tion to the dereddened data, although a slightly higher value 
of the reddening (E(B-V)=0.48 mag) had to be adopted. 



5.5.3 Comparison with the Analysis of Leuenhagen & 
Hamann 1998 

Leuenhagen & Hamann (1998) present fits to their 1 A- 
resolution optical spectrum that are not dissimilar in quality 
from our own, although they fit a smaller spectral range. The 
effective temperature of 40 kK derived with the ISA- Wind 
modelling code, is higher than their 35 kK. Our mass-loss 
rate is similar to the one derived by Leuenhagen & Hamann 
(once the different distances are taken into account), while 
our final wind velocity is more than twice as large. Leuen- 
hagen & Hamann (1998) have not determined v^o from the 
C IV line, limiting themselves to fitting the widths of the 
optical lines. We have already seen how these lines are pe- 
culiarly narrow and in apparent disagreement with the value 
indicated by the P-Cygni absorption of C iv A1550. Inspect- 
ing their Fig. 3 we notice that their model optical lines are 
too narrow and suffer from the same problem as our fits with 
a (3=1 velocity law. 

The hydrogen abundance of 10% (by mass) determined 
by Leuenhagen & Hamann (1998) as an upper limit, when 
used in our model, produces a synthetic Ha-He ll line blend 
which is not entirely in disagreement with our observations. 
Therefore, based solely on this fit, we confirm their deter- 
mination. On the other hand, as already noted in Sec. 5.1.1 
the blend's shape might be due to either nebular emission 
or the instrumental profile, since nebular forbidden lines ex- 
hibit the same pedestal. Our ISA- Wind carbon abundance 
is similar to that determined by Leuenhagen & Hamann 
(1998), while the carbon abundance obtained with CMFGEN 
is lower. Our oxygen abundance, obtained with either of the 
two codes, is higher than the abundance derived by Leuen- 
hagen & Hamann (1998). 
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Figure 6. Top panel: the de-reddened ult&avioleli (lUE), opticaf'^.UCLES) and near-infrared (SOFI) spectra of SwSt 1 (dotted line) 
compe*'J<^\i 



SWS) enei^y ^stribution of SwSt^ (sotted line) compare^ tigthe flux distributiogfrgn the best-fitting isg-\^nd stellar model (s^idQ 
line, logA<4.4) and 1200 K and 230 K blackbody curves (solid line, lQgA>4-^X; Ground-based infrared photometry (Allen & Glass 1974; 
Cohen & Barlow 1974) and colour-corrected IRAS fluxes are also pi are symbols). Wavelength units are A and flux units are 
ergs cm~^ A~^. _ ^ 



6 NEBULAR ABUNDANCE ANALYSIS 
6.1 Nebular Line Fluxes 

The nebular spectrum of SwSt 1 is much stronger and 
slightly higher in excitation than that of the PN with [WCIO] 
central stars analyzed by DBS97 and De Marco & Crowther 
(1998). Many nebular He I lines could be measured: amongst 
them we selected those at 4471.5, 5875.7 and 6678.1 A, dis- 
regarding the line at 7065.2 A due to its being dominated 
by collisions and the possibility that it is affected by telluric 
absorption. The [O ll] line complex at 7325 A is made up of 4 
lines with wavelengths 7318.8, 7319.9, 7329.6, and 7330.6 A, 
where the individual components are partly blended. Un- 
fortunately, the blend at 7319 A was affected by an echelle 
order join which removed a part of the flux. In order to ob- 
tain a measurement of the flux in those lines we therefore 
snipped the section containing the join. The snipped line, 
despite missing the central part, still had clearly identiflable 



wings and could therefore be successfully fitted with a Gaus- 
sian. We nonetheless decided to assign a 50% uncertainty to 
it and did not use it in any of the parameter determina- 
tions. The [O ll] lines at 3726.0 and 3728.8 A are resolved 
and completely un-blended with any stellar feature. [O ill] 
lines appear at 4363.2, 4958.9 and 5006.8 A. 

AU four [S ll] nebular lines included in our wavelength 
range could be measured, namely the lines at 4068.6, 4076.3, 
6716.5 and 6730.8 A. The line at 4068.6 A was found to lie 
on an order join. However since abundance results obtained 
with it or with the weaker lines at 4076.3 A were very sim- 
ilar, we decided to treat it normally and only assign it a 
marginally higher uncertainty. To determine the abundance 
of S'^'*', we used A9531.0. The telluric water vapour line at 
9069.126 A affects the measurement of the [S ill] line at 
9068.9 A since its measured ratio to the A9531 line is 0.28, 
lower than the theoretical ratio of 0.403. 

Three lines of singly ionized nitrogen reside in our spec- 
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Table 8. Observed (Fa) and de-reddened (Ia; E(B-V)=0.46) 
nebular line intensities for SwSt 1 from the 1993 UCLES spec- 
trum. The de-reddened H/3 flux is 1.50xl0~^° erg cm~^ . 
Errors are gauged at 10%. The lines' rest wavelengths are from 
Kaufman and Sugar (1986). The spectral resolution is indicated 
at the beginning of each Section. 
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tral range, but as the [N 11] ratio A6583/A6548 was equal to 
1.4, instead of 2.90, we decided that the strong A6583 line 
was saturated (as was Ha) and did not include it in any of 
our determinations. Together with the A6548 line we used 
[N II] A5755. 

To derive the nebular C/H ratio we used lUE obser- 
vations of the semi-forbidden C 11] A2326 and C ill] A1909 
lines. We assume the stellar contribution to A2326 to be 
negligible, while we could de-blend the stellar and nebular 
contributions to C ill AA1906,1909 using the high resolution 
lUE spectrum (Fig. 0). 

Besides the lines listed in Table |^, we identify forbidden 
lines of Fe+ (e.g. [Fe 11] A4606), Fe^+ (e.g. [Fe lll] A4881) and 
possibly Fe^+ (e.g. [Fe iv] A5677). We also identify [Ne ill] 
A3868, [Ar lll] A5191.8 A7135.8 (affected by an order join) 




1902 1906 1910 1914 
Wavelength (A) 

Figure 7. The high resolution lUE spectrum of SwSt 1 showing 
the semi-forbidden stellar-nebular blend of C in] A1906,1909. 



Table 9. The nebular electron temperature and density adopted 
for SwSt 1, along with values from the literature. 

Te (K) log(Ne) E(B-V) Ref. 

10 500 ± 500 4.5 ± 0.2 0.46 this work 

11400 ± 500 5.0 ± 0.1 0.41 dFPV87 

8000 5.0 ± 0.1 0.34 FGC84 



and A7751.1 and [CI ill] A5517.7 and A5537.9. A line list is 
presented in Appendix B. 



6.2 Nebular Temperature, Density and 
Abundances 

In Fig. ^ we plot a nebular diagnostic diagram for SwSt 1. 
From a weighted mean of the diagnostic lines we adopt val- 
ues for the electron temperature and density of 10 500 K and 
31600 cm"^ (log(Ne)=4.5), respectively (Table ||). 

The temperature implied by the [O ill], [S ill] and 
[Ar III] diagnostic ratios is between 11000 K and 10000 K 
for electron densities between 10'^ and 10^ cm~^. The C ill] 
A1909/A1906 ratio of 0.34 indicates logA^e^S.l for Te in the 
range 9000-14000 K. However this line ratio is severely af- 
fected by the blend with the stellar component. De-blending 
by fitting nebular and stellar components to the doublet re- 
sulted in a large uncertainty. 

The diagnostic ratios from singly ionized elements seem 
to agree on a lower density (4.0<log(Ne)<4.5), if we make 
an exception for the [O 11] ratio A3727/A7330. This ratio 
uses lines from two different grating settings and therefore 
might be affected by errors in the flux calibration. How- 
ever if this were the cause of the discrepancy we would also 
have to be wary of the [S ill] A9531/A6312 ratio (which we 
have used above) and the [S 11] A4072/A6725 ratios, although 
these employ the green and red settings and not the blue 
and red settings used by the [O 11] A3727/A7330 ratio. The 
ratios [O II] A3727/3729 and [S 11] A6717/6731 are unreli- 
able because they are at their high density limit. Therefore 
only [S 11] A4072/A6725 yields a reliable density diagnos- 
tic. From [S ill] A6312/A9531, [O iii] A4959/A4363, [Ar iii] 
A7751/A5192, [N 11] A5755/A6548 and [S 11] A4072/A6725 we 
derive Te=10 500 K and log(Ne/cm-^)=4.5. 

Although the HST images point to a somewhat strati- 
fied PN (with the [O ill] image being smaller than the H/3 
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Table 10. Nebular ionic abundances for SwSt 1. 
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He+/H+ 
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0+/H+ 
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0.044 ± 0.004 
1.21(-4)'' ± 0.60{-4) 



1.19(-4) ± 0.24(-4) 

1.99(-5) ± 0.20(-5) 

2.44(-4) ± 0.24(-4) 

1.05(-5) ± 0.10(-5) 

9.77(-5) ± 1.19(-5) 

3.94(-6) ± 0.80(-6) 

2.96(-7) ± 0.30(-7) 

4.23(-6) ± 0.80(-6) 

=1.21x10-* 

from Aitken and Roche 
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Figure 8. Diagnostic diagram, log(Ne) vs Te, for SwSt 1. The 
solid lines are obtained from ratios of ions in the first stage of 
ionization, the dashed lines are for ions in the second ionization 
stage. 



one) we feel inclined to average the results obtained from 
the singly and doubly ionized elements. SwSt 1 is a young 
and compact PN and we therefore feel unjustified in using 
different electron density and temperature combinations for 
the singly and doubly ionized regions, a choice that is more 
appropriate for more extended, higher ionization objects. 

Flower et al. (1984) also derived log(Ne)=5.0 for SwSt 1 
from the high resolution lUE spectrum of the C ill] 
AA1907,1909 diagnostic lines. They also quoted an upper 
limit for the electron density of 2x10* cm~"^, derived from 
the [O ll] A3729/A3726 ratio (which is almost identical to 
ours), however they argued that this ratio is too close to 
the high density limit to be reliable. Their temperature de- 
termination, using the [O ill] A5007/A4363 ratio, was 8800 
K (for E(B-V)=0.70), not inconsistent with our own value 
(although we employ a lower reddening, the effect is not ma- 
jor over such a short baseline). However their determination 
using the [O ill] A5007/A1663 ratio implies a higher elec- 
tron temperature (15 000 K). If we ascribe the discrepancy 
to the use of two different spectra in the second determina- 
tion (lUE and optical), which might have imperfect absolute 
flux calibrations, we can reconcile our own determinations 
of the temperature and density for the doubly ionised ele- 
ment with theirs, de Freitas Pacheco and Veliz (1987) de- 
rived log(Ne/cm"3)=5.0 and Te = (11400 ± 500) K, from 
[S II] AA4068,76/AA6716,31, [O iii] A4363/(A4959-HA5007) 
and [N ll] A5755/AA6548,84. Their diagnostic diagram is not 
inconsistent with our own. 

Our derived abundances are listed in Tables ^ and |ll| 
where they are compared with the mean values for Type-I 
and non Type-I PN from Kingsburgh & Barlow (1994). The 
quoted errors derive from errors on the fiux measurements 
and do not take into account errors on the electron temper- 
ature and density. The C/H number ratio of 2.4x10"* is 
not as high as for CPD-56°8032 (63x10"*; DBS97), He 2- 
113 (50x10-*; DBS97) or M 4-18 (12x10"*; De Marco & 
Crowther 1998) and lower even than the mean PN value as 
determined by Kingsburgh and Barlow (1994; 5.49x10"*). 
Additionally, with C/O=0.94, SwSt 1 qualifies as an 0-rich 
PN Flower et al. (1984) obtained C/O=0.72, while de Fre- 
itas Pacheco and Vefiz (1987) derived C/O=0.54. The C/O 



ratio is consistent with the presence of a silicate emission 
feature at 10 /xm (see Aitken & Roche 1982 and the SWS 
spectrum plotted in Fig. ^). 

To obtain the N/H number ratio we multiplied the 
N"'"/H"'" ratio by an ionization correction factor of 1.05, ob- 
tained from the ratio of total oxygen to singly ionized oxy- 
gen abundances, in accordance with Kingsburgh and Barlow 
(1994). 



7 HST IMAGES 

Perek and Kohoutek (1967) quote an optical angular diam- 
eter of 5 arcsec, while Aaquist and Kwok (1990) measured a 
5 GHz radio diameter of only 1.3 arcsec. No resolved optical 
image of SwSt 1 was found in the literature. SwSt 1 was 
observed in 1992 by HST in the H/3 and [O ill] A5007 lines, 
as part of a WF/PC snap-shot survey (two exposures per 
filter). These four pre-COSTAR images were already cali- 
brated by the pipeline process (details can be found in the 
HST data handbook). DBS97 used the stellar [O ill] images 
as a point spread function to de-convolve their images of 
CPD-56°8032 and He 2-113. However the nebula around 
SwSt 1 is more highly ionised than the ones around CPD- 
56°8032 and He 2-113, as revealed by the prominent [O ill] 
lines observed in our spectrum. Hence the [O ill] image could 
not be used as a point spread function to de-convolve the H/3 
image. We therefore used the program Tiny Tim (V4.o|^ to 
create theoretical PSFs for the H/3 and [O ill] images (the 
theoretical H/3 PSF is shown in Fig. |^ together with the raw 
images) . 

The de-convolved images of SwSt 1 (Fig. 0) reveal a 
broken ring of diameter ((1.3x0.9) ± 0.2) arcsec, elongated 
in the east-west direction, in agreement with the radio im- 
ages of Kwok et al. (1981) and Aaquist and Kwok (1990), 
but smaller than the 5 arcsec quoted by Perek and Kohoutek 
(1967). Although the images are of poor quality, this struc- 
ture could be a ring whose plane is at about 30° from the 
line of sight. If we take 1.3 arcsec as the true diameter, a 
distance of 2.0 kpc and the expansion velocity of 21 km s"'^ 
would imply a dynamical time-scale of 290 yr for the outer 



^ Ti ny Tim is supported bv the Space Teles cope Science Insti- 
tute; http: / /www.stsci.edu/software/tinytim/ . 
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(c) 

* 




Figure 9. HST images of SwSt 1; a) The average of ttie two raw H/3 images; b) the average of the two raw [O III] images; c) the synthetic 
point spread function generated by the program Tiny Tim; the FOV is 6 arcsec, north is towards the top, east to the left. 



Figure 10. De-convolved H/3 (a) and [O III] (b) HST images of SwSt 1; The images are displayed on a log scale; the FOV is 2.6 arcsec, 
north is towards the top, east to the left. 



edge of the ionized material (although we do not know how 
far the neutral shell extends). 

The [O III] image appears smaller than the H/9 image 
(about 75% of the area), as can be seen from overlaying 
the normalised contours of the [O ill] image onto the grey- 
scale plot of the H P image (Fig. This is consistent with 
ionization stratification of the PN. Further confirmation of 
the stratified structure of this PN comes from analysing the 
nebular line widths as a function of ionic charge (Fig. 
After de-convolving the FWHM of individual lines with the 
respective instrumental profiles, nebular forbidden lines aris- 
ing from singly ionised ions (e.g. [N ll] A6548 and [S ll] 
A6717), as well as lines coming from neutral ions (e.g. [O l] 
A6300) are systematically broader (FWHM~45-60 km s"^), 
while lines arising from doubly ionised ions (e.g. [O ill] A5007 



and [Ar ill] A7750) have widths similar to the Balmer lines 
(FWHM(H/3)=40.2 km s"^) or narrower (FWHM~30-40 
km s"'^; see Fig. This is consistent with the lines aris- 
ing from higher ionization stages being formed in an inner, 
slower part of the PN, and the lines from lower ionization 
stages being formed further out in an accelerating part of the 
flow. The Balmer emission, coming from the entire ionised 
region, might have been expected to produce broader lines, 
although it might be more intense in the inner, denser region 
of the PN, yielding narrower profiles. 



8 NEBULAR MODELLING 
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We carried out a Zanstra analysis using blackbody as 
well as Kurucz ATLAS9 (Kurucz 1991) model atmospheres 
with log(g)=4.0, the minimum gravity tabulated for atmo- 
spheres of the required effective temperatures. We derive a 
blackbody H I Zanstra temperature of (30 300 ± 500) K, 
while the ATLAS9 model grid yields an effective tempera- 
ture of (34 000 ± 500) K. The radius predicted for the scaled 
blackbody is 2.7 R0 at the adopted distance of 2.0 kpc, from 
which a luminosity of 3550 L0 can be derived. A similar ex- 
ercise for the Kurucz atmospheres yields a radius of 2.3 R0 
and a stellar luminosity of 6300 Lq. 

de Freitas Pacheco and Veliz (1987) derived a black- 
body H I Zanstra temperature of 32 000 K. A comparison 
by Flower et al. (1987) of the stellar ultraviolet and neb- 
ular free-free radio continuum fluxes yielded a blackbody 
H I Zanstra temperature of 36 000 K, although they found 
that a 30 000 K blackbody better fitted their optical and UV 
de-reddened stellar energy distribution. 



Figure 11. De-convolved H/3 HST image of SwSt 1. Over-plotted 
are the contours from the [O in] image (Fig. |lO[). 



8 
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Figure 12. Comparison of selected nebular lines in the AAT- 
UCLES spectrum of SwSt 1. The lines are background subtracted 
and normalised to the peak intensity of the [Ar III] A7750 line. 



8.1 Zanstra Temperatures 

Following the formula of Milne and AUer (1975), for a 5 GHz 
flux of 243 mjy (the weighted average of Kwok et al. 1981, 
Milne and AUer 1982 and Aaquist and Kwok 1990, see Sec- 
tion ^.3| ), we predict an intrinsic H/3 flux of 8.07x10"^^^ 
ergs cm~^ s~^. This H/3 flux agrees within the uncertainties 
with the fluxes of Flower et al. (1984) for their adopted red- 
dening of c(H/3)=0.68. From these de-reddened H/3 fluxes 
and the distance of 2.0 kpc for SwSt 1 we derive a flux 
of hydrogen-ionizing photons for the star corresponding to 
log(Q,(s-i))=46.90. 



8.2 Modelling Strategy and Results 

Photo-ionization modelling of PN is the best way to test 
the internal consistency of an empirical analysis. The Har- 
rington photo-ionization code applied here (Harrington et 
al. 1982) assumes that the nebula can be represented by a 
hollow spherical shell which is ionised only by a central star, 
and sampled by 60 radial grid points. 

The Kurucz model and blackbody luminosities and the 
distance were fixed at the values derived above. An inner 
nebular angular radius of 0.22 arcsec (approximately 5 pix- 
els) was determined from the HST observations, and scaled 
to the adopted distance (2.0 kpc) to yield an inner radius 
of 0.0086 pc. The clumping and vacuum filling factors were 
varied so as to reproduce the de-reddened H/3 fiux. The elec- 
tron density was kept constant throughout the PN, fixed at 
the empirically derived value, except in the case of the WR 
fiux distribution, for which this value overestimated the H/3 
fiux even for extremely low values of the filling factor, a fact 
that was judged non-physical. In theory, such a high density 
PN is likely to be optically thick, such that a correct mod- 
elling procedure should make sure that the model reaches far 
enough to include the edge of the Stromgren sphere. This 
is achieved by increasing the thickness of the PN shell un- 
til the the PN recombines at the outer grid points. On the 
other hand, following this procedure would mean that the 
thickness of the ionized shell and the corresponding H/3 fiux 
would be much larger than observed, a fact that cannot be 
compensated for by reducing the filling factor within sen- 
sible limits. We therefore used the observed values for the 
inner and outer radii and accepted an optically thin PN. A 
non-constant density profile was not attempted. The efi^orts 
of De Marco & Crwother (1999) on the PN M4-18 did not 
produce better results when a variable density profile was 
adopted. Additionally, for the compact PN SwSt 1 the lack 
of a high signal-to-noise ratio image means that the density 
profile cannot be constrained. 

As a starting point we adopted the empirically derived 
abundances. Lines of oxygen and carbon were then fitted by 
changing their abundances. Carbon and oxygen are impor- 
tant coolants and changing the abundances not only varies 
the line strengths, but also the nebular electron tempera- 
ture. When a compromise was reached on the strength of 
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Table 11. Empirical and modelled nebular abundances for SwSt 1 compared with 
mean values for Type-I and non Type-I PN from Kingsburgh &; Barlow (1994). 



± VOj LIU 








A Vai 1 in *n Ti r'd 
ii U UliLlaylitjt; 




Empirical 


Model*" 


Type-I 


non Type-I 


He+ /H+ 


0.044 ±0.004 


0.04 


0.129 ± 0.037 


0.112 ± 0.015 


log(C/H)+12 


8.38 ± 0.25 


8.40 


8.48 ± 0.30 


8.81 ± 0.30 


\oz{N/H)+12 


7.32" ± 0.10 


7.99 


8.72 ± 0.15 


8.14 ± 0.20 


log(0/H)+12 


8.41 ± 0.20 


8.15 


8.65 ± 0.15 


8.69 ± 0.15 


log(Ne/H)+12 


8.01 ± 0.07 


7.88 


8.09 ± 0.15 


8.10 ± 0.15 


log(S/H)+12 


6.65 ± 0.10 


6.62 


6.91 ± 0.30 


6.91 ± 0.30 


C/O 


0.94 


1.79 


0.68 


1.32 


N/O 


0.08 


0.08 


1.17 


0.28 



"ICF(N)=O/O+=1.05, N/H=IFC(N)xN+/H+, following Kingsburgh & Barlow 1994 
''Abundances derived from the photoionization model, using the ISA-Wind WR 
stellar atmosphere as input (Table lid) . 



these lines and the electron temperature, lines of sulphur, 
nitrogen and neon were fitted by varying their abundances. 
The abundances of other unobserved species were left at the 
mean PN value as determined by Kingsburgh and Barlow 
(1994) for a set of 80 southern PN, or at their solar values. 
When the ionization balance of the element in question is 
very different from the observations it is arbitrary to decide 
which lines should be fitted and which should be ignored. In 
the case of sulphur, we used the strong [S ill] line at 9532 A, 
ignoring [S ll] lines. For oxygen, we fitted a compromise of 
the [O ii] lines AA3726,29 and A7325, forsaking [O iii] lines. 
For neon, the [Ne ll] line at 12.8 /xm was used, while for 
carbon a compromise between the C ll] A2326 and C ill] 
A1908 lines was used. We note however that because none 
of the stellar flux distributions reproduces the observations 
to an acceptable precision, we found it unnecessary to fit 
individual spectral line fiuxes to high accuracy. 

In Table ^ we compare the abundances determined 
from modelling with the empirical values, as well as the pre- 
dicted line fluxes with the de-reddened observations. None of 
the stellar atmospheres used reproduces the observed char- 
acteristics of SwSt I's PN to an acceptable degree of ac- 
curacy. In the case of the WR flux distribution, we decided 
that lowering the electron density to about half the observed 
value was preferable to having an extremely low filling fac- 
tor (~0.02, too low even if the PN is actually a ring as 
the HST image might suggest). With this value of the elec- 
tron density, a more acceptable filling factor of 0.76 was 
needed to reproduce the H/3 flux. However, this implies an 
optically thin PN. The overall ionization balance is over- 
estimated except that of sulphur, which is underestimated. 
Interestingly, a cooler WR stellar atmosphere, also failed to 
reproduce the ionization balance. We tried to use the 35kK 
model which reproduces most of the stellar spectral features 
except the C iv/Clll ionization balance (see Sect. [S.5| ). This 
model severely underestimates the nebular ionization bal- 
ance, while at the same time having a similar problem with 
the II/3 flux and the flUing factor, as outlined above. We 
therefore believe that the failure of the PN model is not due 
to the hot WR atmosphere. 

The blackbody and Kurucz stellar flux distributions un- 
derestimate the ionization balance of carbon, but overesti- 
mate that of oxygen. They give mixed results in the case of 
neon and sulphur. With particular reference to oxygen, we 
point out that the [O ill] lines used are extremely suscep- 



Table 12. Comparison of results from photo-ionization models 
for SwSt 1 with observations (I(H/3)=100), using ISA-Wind non- 
LTE, blackbody and Kurucz plane-parallel energy distributions. 



Parameter 


Observed 


WR 


Blackbody 


Kurucz 




- 


40.0 


30.3 


34.0 


log(L/L0) 


- 


3.94 


3.74 


3.80 


R(R0) 


- 


2.0 


2.7 


2.3 


log Qo(s-i) 


46.90 


47.53 


46.90 


46.90 


log(I(H/3)) 


-10.1 


-10.1 


-10.1 


-10.1 


Ne(cm~3) 


31620 


15000 


31620 


31620 


e 


- 


0.76 


0.15 


0.18 


Rinner (P*--) 


0.0086 


0.0086 


0.0086 


0.0086 


Router (pc) 


0.0252 


0.0252 


0.0252 


0.0252 




- 


0.021 


0.009 


0.015 


r(H I) 


- 


0.69 


1.3 


2.2 


He/H 


0.04 


0.04 


0.064 


0.15 


C/HxlO'» 


2.51 


2.51 


4.68 


2.51 


N/HxlO^ 


1.99 


9.86 


4.11 


2.13 


O/HxlO'' 


2.57 


1.40 


1.90 


0.78 


S/HxlO"^ 


4.17 


4.00 


6.00 


4.60 


Ne/HxlO"' 


9.80 


7.60 


8.80 


7.90 


1909 C III] 


71 


68 


33 


14 


2326 C ii] 


85 


19 


76 


208 


3726 [On] 


30 


57 


39 


38 


3729 [On] 


18 


23 


14 


14 


3868 [Ne iii] 


0.5 


0.0 


4.8 


0.06 


4068 [Sii] 


2.2 


2.5 


2.5 


3.7 


4076 [Sii] 


0.85 


0.27 


0.82 


1.2 


4363 [Oiii] 


0.3 


0.09 


0.4 


0.02 


4471 He I 


2.5 


2.3 


3.3 


1.2 


5007 [O III] 


29 


204 


62 


2.2 


5755 [Nil] 


3.2 


2.2 


1.9 


3.1 


5876 He i 


7.1 


6.7 


9.2 


3.6 


6300 [Oi] 


1.5 


0.0 


0.06 


0.07 


6363 [Oi] 


0.5 


0.0 


0.02 


0.02 


6548 [Nil] 


31 


30 


27 


32 


6678 He I 


1.7 


1.8 


2.5 


1.0 


6717 [Sii] 


0.7 


0.2 


0.8 


1.1 


6731 [Sii] 


1.4 


0.4 


1.9 


2.5 


7325 [On] 


28 


18 


20 


23 


9532 [S III] 


91 


84 


92 


94 


128 000 [Neil] 


68 


67 


65 


67 


C2+/C+ 


0.98 


9.8 


0.23 


0.17 


02+/0+ 


0.043 


0.81 


0.18 


0.088 


S2+/S+ 


14 


9.7 


10.7 


15 


Ne2+/Ne+ 


0.040 




0.063 


0.004 


Te(N+)(K) 


10 500 


10 500 


9300 


10 800 
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tible to small temperature changes, in particular it appears 
that for central star effective temperatures above ~38kK 
the [O III] lines increase dramatically in flux (Miriam Pena, 
priv. comm.). For these stellar flux distributions we were 
able to use the electron density determined from observa- 
tions, although, once again, we had to accept an optically 
thin PN, so as to obtain a sensible value for the filling fac- 
tor. An optically thin PN is however more in disagreement 
with an electron density of 31600 cm~'^. The three stellar 
atmospheres are compared in Fig. 

The nebular mass implied by all the stellar flux distribu- 
tions is very low compared to the canonical value of 0.3 
for optically thin PN. As for the optical depth properties 
of this PN, it is unusual to have an optically thin PN with 
an electron density as high as we encountered. On the other 
hand the fact that the Zanstra analysis returns values of 
the effective temperature lower than determined with an ad 
hoc modelling of the stellar spectrum with a non-LTE code, 
might suggest that the If/9 flux is not representative of the 
hard photon population, i.e. that the PN is indeed optically 
thin. 

In conclusion, no matter the method chosen for the 
modelling effort, we have found no way of fitting the ob- 
served properties of this PN to a satisfactory accuracy. We 
take this to mean that the modelling tool being used is inad- 
equate for the PN of SwSt 1. In particular the simple shell 
geometry might be far from the truth of this young PN. 
Nebular modelling with our CMFGEN models atmosphere as 
input was not attempted. Overall the two fiux distributions 
compare reasonably well, with ISA- Wind harder far UV flux 
being compensated by CMFGEn's higher effective tempera- 
ture. CMFGEN would therefore not succeed in reproducing 
the [Olll] lines. 



9 SWST 1 BETWEEN 1895 AND 1993 

The warm oxygen-rich dust revealed by the 10 silicate 
emission feature, and the high nebular density, indicate that 
SwSt 1 has recently made the transition from the AGB. 
This could indicate that the pulse that caused SwSt 1 to 
become H-deflcient might have occurred in the last century, 
i.e. after it was first observed in 1895. Since variability of the 
stellar and nebular lines has been noted by Carlson & Henize 
(1997) and by de Freitas Pacheco & Veliz (1987), we decided 
to investigate the variability issue. Photographic spectra of 
SwStl were examined by taking a cut through a scan of the 
plate that appeared in the original paper. We then applied a 
rough wavelength calibration. Obviously, some fidelity is lost 
in this process but the spectra are still useful qualitatively. 

SwSt 1 was discovered to be an emission line star by 
James E. Keeler in 1895 (Fleming 1895). He reported Balmer 
lines in emission and compared the spectrum to that of 77 
Car. Although we could not find his spectrum in the lit- 
erature, we examined a scan from another plate taken in 
1896 which shows only 3 emission lines (II/3 and [O ill] 
AA4959,5007) superimposed on a stellar continuum. On the 
basis of this spectrum, SwSt 1 is referred to as P-Cygni 
class by Fleming in her log-book. Cannon (1916), referring to 
these or similar spectra, also assigned SwSt 1 to the P-Cygni 
class. These early spectra of SwSt 1 are very low dispersion 
prism spectra. Although P-Cygni profiles had already been 



recognized in other stars, it is doubtful that they could be 
resolved on these plates. It is more likely that SwSt 1 was 
classified as P-Cygni on the basis of bright Balmer lines on 
top of a stellar continuum. Cannon (1921) mentions the pres- 
ence of [O III] AA 4959,5007 emission in SwSt 1. She notes 
that these lines are associated with gaseous nebulae (even 
though they weren't identifled until 1928). 

Subsequently, SwSt 1 was observed spectroscopically on 
six distinct occasions: by Payne Gaposchkin (1938, source 
spectra unknown); in 1940 by Swings & Struve (1940; in the 
range 3600-6680 A at dispersions of 100 and 50 A/mm); in 
1942 by Swings & Struve (1943; 3600-5300 A); in 1962 by 
Carlson and Henize (1979; in the range 3700-6600 A at a 
dispersion of 177 A/mm blue-ward of 5000 A and 94 A/ mm 
red-ward of 5000 A); by AUer (1977, although the date of 
the spectrum is not indicated; the range of the observations 
is 5000-8500 A); in 1976 by Flower et al. (1984; in the range 
4260-6710 A at a resolution of 7 A); in 1985 by de Freitas 
Pacheco and Feliz (1987; in the range 3600-6930 A at res- 
olutions of 6.5 and 0.5 A); in by Leuenhagen & Hamann 
(1998; in the range 3800-5000 A and 5800-6700 A at a reso- 
lution of 1 A) and by us in 1993 (in the range 3700-5400 A 
at a resolution of 0.15 A and in the range 5400-9200 A at a 
resolution of 0.35 A). 

Two trends since 1940 have been suggested in the liter- 
ature: a weakening of the stellar wind lines and an increasing 
degree of ionization of both the stellar wind and the PN. The 
former trend was pointed out by Carlson and Henize (1979), 
who noted the possible disappearance of several weak stellar 
O III and C III lines between 1940 and 1962, along with the 
appearance of [Fe 11] and [Fe ill] lines. They also reported 
the diminution of P-Cygni profiles, from a total of 10 lines 
(as counted by Swings & Struve, (1940)) to 4 lines, finally 
concluding that the wind's density and mass-loss rate must 
have decreased between the two observations. 

We looked at the plate reproductions in the original 
1940 and 1962 papers and we also took a cut along the spec- 
tral direction of scanned-in versions. To within the precision 
allowed by our method, we notice no obvious change. It is 
possible that weak P-Cygni profiles were apparent to Swings 
& Struve in the original plate material and that these fea- 
tures had become weaker by 1962 when observed by Carlson 
& Henize. One might in fact expect the star to have gradu- 
ally lost its hydrogen during its early post- AGB phase. We 
have no reason to doubt the claim by Swings and Struve that 
H7 and He had P-Cygni profiles but not Ha, H/3 or H5. How- 
ever, we would like to point out that Ha is generally the first 
line to develop an emission component and a P-Cygni profile 
in the presence of a wind. Looking at the 1940 spectrum we 
do notice however that the stellar continuum is extremely 
irregular, while in 1960 it is much smoother. Such irregu- 
larity, combined with a variety of emission lines could have 
given the impression of troughs, which could in turn have 
been interpreted as P-Cygni. 

The claim of Carlson & Henize (1979) that [Fe 11] and 
[Fe III] emission appeared sometime between 1940 and 1962 
is also somewhat suspect, since [Fe 11] emission was already 
reported previous to 1940, by Payne Gaposchkin in 1938. 
However by far the most misleading statement about vari- 
ability reported in the literature is the claim of Carlson and 
Henize (1979) that the C ill A5696 and C iv AA5801,12 lines 
are absent from the 1962 spectrum of SwSt 1. These lines are 
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Figure 13. A comparison of the ISA- Wind line-blanketed non-LTE model atmosphere (dotted line) with the 30 300 K blackbody flux 
distribution (dashed line), and a 34 000 K, log(g)=4.0 plane-parallel Kurucz model atmosphere (thick solid line). The H I and He I and 
other important ionization edges are marked. 



prominent in the spectrum of Swings & Struve (1942) and 
also in our own spectrum. The Carlson & Henize (1979) 1962 
spectrum longward of 5000 A comes from a very weakly ex- 
posed plate compared to their other material. On that plate 
only the strong nebular [N ll] A6583 line and Ha are visible, 
with He I A5876 described as being barely visible. The C ill 
A5696 and C iv AA5801,12 lines are several times weaker 
than A5876 in our spectrum and would be well below their 
plate limit. On this basis, the claim that C ill A5696 and 
C IV AA5801,12 are not present in 1962 can be dismissed. 

The only spectral change which seems real, comes from 
the spectroscopy of AUer (1977). From his Figure 4 we see 
that C IV A5470 is much stronger and broader than the 
neighbouring He ll/[Fe ill] blend at 5411 A. In the spectrum 
presented here, the strength of the two lines is comparable. 
However the C iv lines at 5801 and 5812 A are not stronger 
in the AUer spectrum, so that no consistent trend is noticed. 

de Freitas Pacheco & Veliz (1987) advocate a second 
trend, increasing ionization of both the stellar wind and the 
PN. They were not able to see any nebular [Fe ll] lines in 
their spectra despite the far superior resolution of their data 
to earlier spectra, and noticed instead weak lines which they 
attributed to [Fe iv]. Although de Freitas Pacheco and Veliz 
do not show any of the spectral regions containing [Fe ll] 
lines, from the quality of their published spectroscopy we 
suggest that they would not have been able to observe the 
weak [Fe ll] lines detectable in our 1993 spectrum. However, 
since the 1938 and 1962 spectroscopy is likely to have been 
of poorer S/N than their 1985 data, we conclude that [Fe ll] 
lines have indeed weakened and that the fact that we do 
observe them in the 1993 spectrum is only because of our 
higher resolution. Based on the weakening of the nebular 
[Fe ll] lines, we conclude that the PN might have increased 
its level of ionization, although we attach a large uncertainty 
to this statement. However, the argument that the stellar 
wind has increased in ionization is based on the erroneous 



claim by Carlson and Henize (1979) of the absence of the 
C IV lines at 5801 and 5812 A in 1962. 

In conclusion the arguments for a significant evolution 
of the spectrum of SwSt 1 in the last century are weak. The 
PN might have increased its level of ionization slightly, and 
there might have been some variability in the stellar wind 
spectrum, but nothing major has happened to this central 
star in the last 60 years and possibly in the 110 years since 
it was first observed. 



10 DISCUSSION 

From the properties derived for SwSt 1 we can state the fol- 
lowing. Although the possible distance range of 2-4 kpc is 
large, even at a distance of 2 kpc the stellar luminosity im- 
plied by the models points to a central star mass at the high 
end of the PN distribution. As we have seen, the distance 
issue poses a significant obstacle and does not allow us to 
determine the the stellar and nebular parameters with suf- 
ficient accuracy to establish relationships between difi'erent 
objects. Better parallaxes, or expansion distances are the 
only way to get around the problem. Another way would be 
to find [WCL] stars in the Magellanic Clouds and use those 
to characterize the sample. A wide field survey would be 
needed, sensitive enough to detect emission line stars with 
equivalent widths of the order of a few Angstroms and for 8 
m-class telescope follow up spectroscopy. 

The stellar analysis indicates a hotter star than pre- 
viously implied (Leuenhagen & Hamann 1998), and more 
in conformity with the [WC9] spectral classification. The 
weakness of SwSt I's stellar lines can be attributed to a 
relatively low mass-loss rate, although the derived value of 
1.90 X 10~^ Mq yr~^ should be quoted with care because of 
its distance dependence. Also noteworthy is the stellar wind 
oxygen abundance. This determination is reliant on a single 
line of O ill, but no other lines exist in the available spectral 
range which contradict this determination. The new stellar 
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evolutionary calculations of Herwig (2000) show that the 
oxygen abundance of models which account for convective 
overshoot (which begin to account for the severe hydrogen 
deficiency of WR central stars) should have a 10-20% oxygen 
mass fraction. This is in agreement with our oxygen deter- 
mination. Additionally the C/0 mass ratio is in the range 
2-4 (depending on what modelling effort we rely on), show- 
ing third dredge-up enrichment. This is not dissimilar from 
the values obtained for the [WCIO] stars CPD-56°8032 and 
He 2-113 of 3.7 and 4.1, respectively. 

The PN of SwSt 1 is compact and of high density. It 
is likely to be stratified as indicated by the HST H/3 im- 
age being larger than the [O ill] image and by the lines 
from higher ionization stages being narrower. The nitrogen 
abundance is not particularly high, which is surprising if the 
central star is as massive as is implied by other parameters. 
Its other nebular abundances are not particularly unusual 
when compared to the mean values obtained by Kingsburgh 
and Barlow (1994) for non-Type I PN. We confirm again 
that the C/O ratio is smaller than unity (0.7 by mass), indi- 
cating an 0-rich PN. This is interesting in view of the stellar 
C/O ratio of 2-4, indicating active dredge up between when 
the nebular material was ejected the present time. A study 
of the C/O variation as a function of distance from the star 
is underway for SwSt 1, CPD-56°5682 and He 2-113 (De 
Marco et al. in preparation). 

The PN neutral envelope presents a fair amount of 
atomic (e.g. sodium; Dinerstein et al. 1995) and only weak 
molecular H2 emission is observed in our SOFI spectrum. 
We find that modelling this PN with a simple shell geom- 
etry and a smooth electron density fails no matter what 
stellar fiux distribution is used. For a distance of 2 kpc, our 
best fitting model has a luminosity of 8900 , significantly 
higher than the re-radiated IR dust luminosity of 1450 Lq, 
implying that either the dust shell does not absorb and re- 
radiate the stellar luminosity, or else that our stellar model 
is not valid. On this line we also note that the relatively low 
N/H PN abundance ratio does not point to a particularly 
massive central star. Additionally, if wind strength corre- 
lates with luminosity, as is the case for radiationally driven 
winds, the weakness of the wind would also tend to point to 
a lower mass. 

After carefully assessing the spectral variability re- 
ported in the literature by going back to the original publi- 
cations and spectral plates, we cannot fully corroborate any 
of the changes reported in the spectrum of SwSt 1. From this 
we conclude that the star has not changed significantly over 
the past 100 years. Certainly we have not caught SwSt 1 
before the pulse that transformed it into a hydrogen defi- 
cient star. The dynamical time-scale of 150-330 yr refers to 
the edge of the ionized PN, although for such a high density 
object, we might expect this to be only a lower limit, since 
the PN material might not be fully ionized. 
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APPENDIX A: THE NEUTRAL ENVELOPE OF 
SWST 1 

The Na I D lines in the spectrum of SwSt 1 shown in Fig. 1 
indicate a complex neutral envelope structure. Below we 
study this complex line system and compare our findings 
to those of Dinerstein et al. (1995). 

The Na I D line negative LSR radial velocity troughs 
in the spectrum of SwSt 1 are unlikely to be of interstellar 
origin since the Galactic radial velocity curve in the direc- 
tion of SwSt 1 is positive (Fig. Ic) . It is therefore likely that 
these components arise from a neutral envelope enclosing 
the ionized PN. As such they can be a good tracer of veloc- 
ity components. We identified two components per line, with 
radial velocities at (-78 ± 5) km s"'^ and (-38 ± 2) km s~^ 
(Fig. lb). If corrected for the heliocentric radial velocity of 
the nebula, the Na I D line emission components lie at (0 ± 
4) km s~^, while the two troughs are measured to be at (-70 
± 5) and (-28 ± 2) km s"^ The nebular expansion veloc- 
ity determined in Section 3.2 is (20 ± 1) km s"'^. Hence we 



would expect absorption components of the Na D lines asso- 
ciated with the ionized nebula to be found at a heliocentric 
radial velocity of about -20 km s~^. The trough associated 
with the first component is at about -28 km s~^, 8 km s~^ 
more negative than expected, although this discrepancy can 
be almost justified by the relative uncertainties. The com- 
ponent at -70 km s~^ then has to be explained by assuming 
that SwSt I's PN has a neutral velocity component mov- 
ing at a higher velocity than the ionised 'shell' measured by 
the Balmer lines. For the [WCIO] stars CPD-56°8032 and 
He 2-113 (DBS97) there is a correspondence between the 
high velocity components of the Na I D line troughs and the 
shape of the emission lines of [O l] . The same comparison in 
SwSt 1 reveals that fitting the [O l] lines with two Gaussians 
(of which one is forced to have a FWHM=40 km s^^) cannot 
encompass components with FWHM larger than 60 km s~^. 
This could explain a Na I D line component at about -30 
km s~^ (like the observed component at -28 km s~^) but 
not the trough observed at -70 km s~^. In conclusion, the 
absorbing sodium must be associated with a different part 
of the nebula from the coUisionally excited neutral oxygen. 

Dinerstein et al. (1995) also found evidence for circum- 
stellar sodium around SwSt 1. They measured the emission 
component of the Na D lines to be at -7 km s~^, while they 
determined the absorption component to be at -27 km s~^ 
(where their measurements were corrected for a heliocentric 
nebular radial velocity of -20 km s~^). These correspond 
within the uncertainties to our emission component at 
km s~^ and absorption component at -28 km s~^. They 
did not find a double absorption as we do, and indeed their 
Fig. 2 is significantly different from our Fig. ^ (a) , although 
their claimed resolution (0.20 A) is similar to our own. It is 
not possible at present to resolve this discrepancy. However 
they did mention the possibility of the existence of higher 
velocity components in the PN envelope and they compared 
this effect to the fast low ionization outflows (FLIERs) de- 
scribed by Bahck et al. (1993) and Bahck (1994). 
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APPENDIX B: THE OPTICAL AND NEAR-IR 
SPECTRA OF SWST 1 



Table Bl. Lines identified in tlie spectrum of SwSt 1 between 
3700 and 9700 A. 



In this Appendix we present the UCL Echelle spectrum of 
SwSt 1 obtained at the Anglo-Austrahan Telescope in May 
1993 (Figs, 
trum (Fig. 



psD as well as our near-IR NTT SOFI spec- 
details of the data acquisition and reduction 
are presented in Section Due to the strength of the nebu- 
lar lines compared to the stellar lines, we present the rectified 
optical spectrum on a log scale. 

In Table |^ we present a list of the lines we identified in 
the optical spectrum of SwSt 1. In Column 4 we indicate the 
origin of the line (s for stellar, n for nebular, is for interstel- 
lar) , while in Column 5 we list the initials of the authors who 
have identified those lines in the past. A discussion of the 
stellar and nebular compone nts of SwSt I's optical spectrum 
has been presented in Sees. 



5.1 



and 3.1 



The near-IR SOFI spectrum is dominated by strong 
nebular emission lines of Hi (e.g. P/3, Br7) and He I 
(1.083um, 2.058Atm), plus [Slll] 9530A. As reported in 
Sec. O, we have detected H2 (at 2.1218/im {v=l-Q, S(l)) 



for the first time from this young FN. 

This paper has been produced using the Royal Astronomical 
Society /Blackwell Science M^^X style file. 
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Figure B6. The SOFI spectrum of SwSt 1. 
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Table B2. Principal near-IR line fluxes, F^, obtained from NTT- 
SOFI observations of SwSt 1, together with de-reddened fluxes 
relative to I(H/3)=1.50 xlO"'^'' ergcm~^s~^, assuming E(B- 
V)=0.46 mag. Uncertain values are given in parenthesis (due to 
poor atmospheric transmission). 
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This figure "swst_beta_conv_log_lab.gif" is available in "gif" format from: 
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